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Levy walks and scaling in quenched disordered media 
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We study Levy walks in quenched disordered one-dimensional media, with scatterers spaced ac- 
cording to a long-tailed distribution. By analyzing the scaling relations for the random-walk prob- 
ability and for the resistivity in the equivalent electric problem, we obtain the asymptotic behavior 
of the mean square displacement as a function of the exponent characterizing the scatterers distri- 
bution. We demonstrate that in quenched media different average procedures can display different 
asymptotic behavior. In particular, we estimate the moments of the displacement averaged over 
processes starting from scattering sites. Our results are compared with numerical simulations, with 
excellent agreement. 
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Random walks in quenched random environments oc- 
cur in many fields of statistical and condensed matter 
physics [l[ , as they represent the simplest model of diffu- 
sion phenomena and non-deterministic motion. Disorder 
and geometrical confinement are known to strongly influ- 
ence transport properties. In particular, in highly spatial 
inhomogeneous media, the diffusion process is often char- 
acterized by large distance diffusion events, which play 
a crucial role in transport phenomena and can strongly 
enhance them Q. Molecular diffusion at low pressure in 
porous media is dominated by collision withpore walls, 
with ballistic motion inside the large pores Q, diffusion 
in chemical space over polymer chains can be described 
by a distribution of step length with power law behavior 
[J]. In addition, recent experiments on new disordered 
optical materials paved the way to a tuned engineering 
of Levy-like distributed step lengths |5| . These and many 
other processes can often be successfully analyzed using 
the Levy walks formalism Q: The random walker can 
perform long steps, whose distribution is characterized 
by a power law behavior A(r) ~ l/r Q+1 , with a > 0, for 
large distance displacements r. 

An important feature of these experimental settings is 
that the random walk is in general correlated, and the 
correlation is induced by the topology of the quenched 
medium. Diffusing agents moving in highly inhomoge- 
neous regions, where they just experienced a long dis- 
tance jump without being scattered, have a high proba- 
bility of being backscattered at the subsequent step un- 
dergoing a jump of similar size, and this leads to a cor- 
relation in step lengths. While the effect of annealed 
disorder on transport properties in Levy walks is quite 
well understood [7|], the role of correlations in Levy-like 
motion is still an open problem. 

If the motion occurs in low dimensional samples, spa- 
tial correlations in jump probabilities can deeply influ- 
ence the diffusion properties. This was first evidenced in 
models of Levy flights, [8] and more recently discussed in 
one-dimensional models for Levy walks on quenched and 



correlated random environments. The recent studies fo- 
cused, respectively on the mean square displacement in a 
Levy-Lorentz gas [9j , and on the conductivity and trans- 
mission through a chain of barriers with Levy-distributed 
spacings [l(|. Both studies evidenced the effect of cor- 
relations on large scale quantities: the mean square dis- 
placement and the conductivity indeed feature a modi- 
fied scaling behavior with respect to the annealed case, 
and similar results have been obtained for a deterministic 
fractal characterized by holes distributed according to a 
Levy statistics (Tlj . 

We consider here an analogous model. In particular 
we introduce a one-dimensional structure where we place 
scatterers spaced according to a Levy-type distribution 
(Figure [1]) so that the probability density for two consec- 
utive scatterer to be at distance r is 
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where a > and ro is a cutoff fixing the scale- leng th 
of the system. A continuous time random walk |l3[ is 
naturally defined on this structure, i.e. a walker moves 
ballistically (at constant velocity v) until it reaches one of 
the scatterers, and then it is transmitted or reflected with 
probability 1/2. Figure Q] also describes the equivalent 
electrical model, where the resistance between two points 
is defined as the number of scatterers between them [l(| ■ 
A subtle effect of quenched disorder is the dependence 
of the observables on the choice of the walker starting 
site or, equivalently, of the position of the electric con- 
tacts. In particular, on inhomogeneous structures, aver- 
ages over different starting points can provide different 
results with respect to the corresponding local quanti- 
ties [ll], EH- In presence of long tails, local quantities 
are strongly fluctuating and averages are in some sense 
necessary for a global characterization of the structure. 
In our system, two different average procedures can be 
introduced and interestingly, due to Levy statistics, they 
provide different results [9]. 
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FIG. 1: The Levy walk model and the equivalent electric 
problem. Scatterers are placed at positions o, n, rz . . . . The 



spacmgs rj + \ 



between scatterers are Levy distributed 



according to the probability density JT]). 



The first possibility is to consider as a starting site at 
t = any point of the structure and then average over all 
stochastic processes and realizations of the disorder. In 
this situation the behaviour is quite clear @. For a < 1 
the motion is always ballistic, while for a > 1 it can be 
shown that the probability of reaching a site at distance 
r at the first step decays as l/r a (i.e. much more slowly 
than A(r)) [|; hence, the motion can be dominated by 
the first jump (called ballistic peak). In particular, for 
1 < a < 2 the first step provides the major contribution 
to the mean square displacement, i.e. (r 2 (t)) ~ Ct 3 ~ a . 

In experimental realizations however, processes begin 
with a scattering event, corresponding to the entrance 
of particles or light in the system Q and, therefore, av- 
erages should be performed considering only scatterers 
as starting points. In this case the static behaviour of 
the resistance has been calculated in jlOj, while the sit- 
uation in the study of the dynamics of random walks is 
much more obscure. Indeed, the probability of reach- 
ing a distance r at the first step is the same as that 
of any other scattering event. Therefore, the first-jump 
does not determine the behaviour of the means square 
displacement and only a lower bound can be proved [9|: 



(r 2 (t)) > v 2 t 2 P+{t) 



where P + (t) is the proba- 



bility for the first step to be longer than t. An insight 
on the issue would represent a relevant result both for 
a theoretical understanding of quenched Levy processes 
both for guiding future experimental investigations. 

In this paper, using the analytic result for the resis- 
tance lldl . the mapping between resistance and random 
walks 14], a scaling hypothesis for the probability dis- 
tribution P(r,t) [15| and the reasonable hypothesis of 
"single long jump" , we derive an analytic expression for 
the asymptotic behaviors of the characteristic length l(t) 
of P(r,t) and for the mean square displacement (r 2 (t)), 
when averaged over the scattering points. The predicted 
behaviors are 



£(t) 



f i+<* if < a < 1 
ti if 1 < a 
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(r 2 (t)) 



2+ i+°° if < a < 1 
ii~ Q if 1 < a < 3/2 
t if 3/2 < a 
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Our predictions are compared with extensive numerical 
simulations, with excellent agreement. We note that for 
< a < 3/2 the system does not follow the standard 
behaviour (r 2 (t)) - £(t) 2 . 

Let us introduce our argument. In a one dimensional 
system the most general scaling hypothesis for P(r,t) is: 



P(r,t)=r 1 (t)f(r/£(t))+g(r,t) 
with a convergence in probability 



(4) 



lim 

t—*oo 



\P(r,t)-£- 1 (t)f(r/e(t))\dr = (5) 



The integration cut off is provided by the fact that the 
walker in a time t covers at most a distance vt (v is the 
velocity). The leading contribution to P(r,t) is hence 
£~ x (t) f (r I i(t)) which is significantly different from zero 
only for r < £(t). The subleading term g(r,t), with the 
condition lim t _ J . 00 J \g(r, t)\dr = 0, describes the behavior 
at larger distances, i.e. £(t) <C r < t. Notice that, if 
g(r, t) does not vanish rapidly enough, it can nevertheless 
provide important contributions to (r 2 (t)). 

The scaling of the correlation length £{t) is given by: 



£{t) - t d °' 2 



(6) 



where we denoted the exponent as d s /2, in order to 
obtain P(0,t) ~ t~ ds / 2 , i.e. the standard behavior 
of the return probability for random walks on graphs 
[III. In [O, da] it is shown that, in the equivalent elec- 
tric network problem, the resistance R(r) between sites 
at distance r corresponds to the random walk quantity 
lim w _>o P{ r i where P(r,uj) is the time Fourier trans- 
form of P(r,t) -P(0,t). 

Since g(r, t) is relevant only in the regime r 3> £{t) and 
the resistance is evaluated at a finite distance r and for a 
diverging characteristic length £(t) (i.e. uj — > 0), it is clear 
that g{r,i) does not provide significant contributions to 
R(r). Therefore, analogously to [O, [IB], the scaling of 
the resistance is given by 



R(r) 



(7) 



The stationary problem and the calculus of the resistance 
are a much more simple task than the direct dynamical 
study of random walks, so that the asymptotic behavior 
of R at large distances can be calculated analytically 
obtaining R(r) ~ r a for a < 1 and R(r) ~ r for a > 1. 
Plugging this result into equations (J5][TJ) we obtain the 
asymptotic behaviour of £{t) ©• 

The explicit expression of the mean square displace- 
ment reads 

r-vt f-vt 

{r 2 {t))= £-\t)f(r/£(t))r 2 dr+ g(r,t)r 2 dr. (8) 
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Anomalies with respect to the usual behavior (r 2 (t)) ~ 
£ 2 (t), can be present for two reasons. First, the second 
term in (jHJ can be dominant with respect the first one 
(this happen e.g. when averaging over any starting site 
[9( ) . A second more subtle anomaly shows up if the scal- 
ing function f(x) decays too slowly for x — > oo, as in 
the case of the annealed model [7fl. We will show that, 
depending on the value of a, our system presents both 
situations. 

Since the anomalies are determined by the regime 
where r ^> i(t), we expect that the dynamics is domi- 
nated by a single long jump of length r. In averages over 
all starting sites the probability of long jumps is much 
higher at the first step, hence the single long range event 
occurs at t = 0. Here on the other hand, it can happen, 
with equal probability, at any scattering site. In particu- 
lar, neglecting the possibility of multiple "long jumps" 



we obtain, for r > £(t), P(r,t) ~ N(t)/r 



l+a 



« 1, 



where N(t) is the number of scattering sites visited by 
the walker in a time t and l/r 1+a is the probability for a 
scatterer to be followed by a jump of length r. Discarding 
the single long jump, the distance crossed by the walker 
in time t is of order £(t). According to the results in [Toj . 
the number of scattering sites visited in this time is the 
resistance of a segment of length £(t), i.e. £(t) a for a < 1 
and £{t) for a > 1. This implies that N(t) ~ for 
a < 1 and N(t) ~ i 1/2 for a > 1. 

For a < 1 and r 3> £(t), we estimate P(r,t) as 



P(r,t) ~ 



£(t) G(i)) 



-(l+a) 



(9) 



Hence the scaling function f(x) features a long tail for 
large x decaying as x~ 1 ~ a . On the other hand for a > 1 
and r 3> £ (i) we obtain 



P{r,t) ~ t*. 
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-(l+a) 



•g(r,t) 

(10) 



Note that here g(r,t) provides a subleading contribution 
to P(r,t). Indeed lim^oc J^L fl"(r,i) = (a > 1). 

We note that the contribution to (r 2 (t)) of lengths r < 
£(t) is always of order £(t) 2 , while, at larger distances, 
the dominant contribution is provided by probabilities 
(j9l fl0|) . The contributions coming from these tails are, 
for a < 1 




FIG. 2: (Color on-line) Mean-square displacements as a func- 
tion of time for different values of a, i.e. symbols. Theoreti- 
cal predictions given by equations ((3J are represented by thin 
(red) lines. Theoretical values of the exponent as a function 
of a is compared with the fitted values in the inset. 




FIG. 3: (Color on-line) Montecarlo evaluation of the probabil- 
ity density rescaled by £(t) as a function of r/t(t) for a = 0.3. 
The dashed black line evidence that for r > l(t) P(r, i) fea- 
tures the behavior described in |[9j- Spikes correspond to 
ballistic peaks at r — t, which provides the subleading contri- 
bution (r 2 ) ~ t 2 ~ a to the mean square displacement. 
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and for a > 1 



.l+a 



dr ~ f 



dr ~ t 2 ' 



(11) 



(12) 



The contribution (fTTj) for large times is always greater 
than £ 2 (t) , while (fl~2"|) is dominant with respect £ 2 (t) 
only for a < 3/2. The overall behavior of the mean 
square displacement is therefore given by ©• 



Analogously we estimate the moments (r p (t)), p > 0, 

if a < 1, p < a 



(r p (t)) 



p(l+°)-° 2 

t ITS 

ti - £(t)P 



if a < 1, p > a 

if a > 1, p < 2a - 1 

if a > 1, p > 2a - 1 



(13) 



We now compare our finding with extensive numeri- 
cal simulations. We set the cutoff for the length scale 
of the system, rp =0.1 and the velocity in the ballistic 
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FIG. 4: (Color on-line) Montecarlo evaluation of the probabil- 
ity density rescaled by l(t) as a function of r/£(t) for a = 1.3. 
The black lines evidence that for r > £(t) P(r, t) features the 
behavior described in (|10[) . 

stretches v = 1. In Figure [2] we plot the mean square 
displacements as a function of time. The numerical re- 
sults obtained with Montecarlo simulations (symbols) are 
compared with the behavior predicted by formula ([3]), 
represented by thin (red) lines. An excellent agreement 
is present in the whole range of a's. 

In Figure [3] we plot the probability density P(r,t) for 
a = 0.3 at different times. The scaling hypothesis (@1 is 
well satisfied if £{t) grows as 

We also verify the 

behaviour for large r predicted by equation (JSJ. Figure 0] 
evidences that, for 1 < a < 3/2, P(r,t) presents at large 



distances a different behaviour. The probability always 
decreases as (r/£(t))~ 1 ~ a but the coefficient depends on 
time, vanishing as it 1 ^ )/ 2 . In conclusion, the numerical 
simulations confirm the single long jump hypothesis that 
we have introduced in the calculation of P(r, t) both for 
< a < 1 and for a > 1. 

In this paper, we have presented an analysis of a one 
dimensional Levy walk in quenched disorder with corre- 
lation induced by the topology. Using the mapping with 
the equivalent electric problem and the scaling hypothe- 
sis for the random walk probability, we have proposed a 
"single jump" approximation leading to a new estimate 
for the mean square displacement and for its moments, 
when averaged over the scattering points, as a function 
of the parameter of the Levy distribution for the spac- 
ings. Our hypothesis is able to capture correctly the slow 
decaying tail effect of the P(r,t), leading to an explicit 
difference between the scaling length £ 2 (t) and (r 2 (t)), 
which is different from the analogous result originating 
from the contribution of the ballistic peak, observed in 
previous works. Interestingly, in experimental settings 
with tuned Levy disorder [j| , the average over the scatter- 
ing points appears to be the relevant quantity to compare 
with experimental results. Our findings open the way to 
an extension of analogous arguments to higher dimen- 
sional disordered samples, where the effects of quenched 
disorder is still unclear, in view of future experiments. 
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J. Bertolotti, S. Lcpri, R. Livi, K. Vynck and D.S. 
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project PRIN 2008 Non linearity and disorder in classical 
and quantum processes. 



[1] J. P. Bouchaud and A. Georges, Phys. Rep. 195 127 
(1990); D. ben-Avraham and S. Havlin, Diffusion and 
Reactions in Fractals and Disordered Systems, Cam- 
bridge University Press, 2004 

[2] R.Klages, G. Radons and I.M. Sokolov (Eds.) Anomalous 
Transport: Foundations and Applications (Wiley, VCH 
Berlin), 2008. 

[3] P. Levitz, Europhys. Lett. 39 6593 (1997) 

[4] I.M.Sokolov, J.Mai and A. Blumen, Phys. Rev. Lett. 79 
857 (1997) M. A. Lomholt, B. van den Broek, S.-M. J. 
Kalisch, G. J. L. Wuite and R. Metzler, PNAS 106 8204 
(2009) 

[5] P. Barthelemy, J. Bertolotti and D.S. Wiersma, Nature 

453 495 (2008). 
[6] M. F. Shlesinger and J. Klafter, Phys. Rev. Lett., 54 2551 

(1985), J. Klafter, A. Blumen, G. Zumofen, and M. F. 

Shlesinger, Physica A, 168 637 (1990). 
[7] T. Geisel, J. Nierwetberg and A. Zacherl Phys. Rev. 

Lett. 54 616 (1985), M. F. Shlesinger, G. M. Zaslavski 

and J. Klafter Nature, 363 31 (1993), G. Zumofen and 

J. Klafter Phys. Rev. E 47, 851 (1993). 



[8] H.C. Fogedby, Phys. Rev. Lett. 73 2517 (1994), R. Kut- 
ner and P. Maass, J. Phys. A: Math. Gen. 31, 2603 
(1998), M. Schulz, Phys. Lett. A , 298, 105 (2002). 

[9] E. Barkai, V. Fleurov, J. Klafter, Phys. Rev. E 61 1164 
(2000). 

[10] C.W.J. Beenakker, C.W. Groth, A.R. Akhmerov, Phys. 

Rev. B 79, 024204 (2009). 
[11] R. Burioni, L. Caniparoli, S. Lepri and A .Vezzani Phys. 

Rev. E 81, 011127 (2010) 
[12] R. Burioni and D. Cassi, J. Phys. A 38, R45-R78 (2005). 
[13] E. W. Montroll and G. H. Weiss, J. Math. Phys. 6 167 

(1965), G. H. Weiss and R. J. Rubin, Adv. Chem. Phys. 

52 363 (1983), J. W. Haus and K. W. Kehr, Phys. Rep. 

150 263 (1987). 
[14] P.G. Doyle and J.L. Snell, Random Walks and Electric 

Networks, The Mathematical Association of America, 

Inc., Washington, D.C. (1999). 
[15] M.E. Cates, J. Physique 46, 1059, (1985). 
[16] S. Alexander and R. Orbach, J. Physique Lett. 43 L62 

(1982). 



